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Interferometic tomography in the presence of an opaque object has been 
investigated. The developed iterative algorithm does not need to augment the 
missing information. It is based on the successive reconstruction of the 
difference field, the difference between the object field to be reconstructed 
and its estimate, only in the defined region. The application of the 
algorithm results in stable convergence. 
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The interest in tomographic reconstruction arises tn various fields including 
interferometric flow visualization,* - ^ medical imaging,^ electronmicroscopy, ** 
and nondestructive testing. ^ The reconstruction in these fields frequently 
encounters ill-posed problems of incomplete projection and limited view angle, 
which pose formidable challenges of restoring distorted images. This is 
especially true in optical fluid-flow measurements, where physical 

constraints (i.e. , test model and test-section enclosure) limit the angular 
scanning or block part of the probing rays. This paper discusses a new 
approach for reconstructing a field from incomplete projections as depicted in 
Fig. 1. Discussion will be developed in terms of interferometric recon— 
truction of continuous flow fields; however, its extension to analogous fields 
is self-evident. 

In interferometry, the projection data which represent fringe order 
numbers are the line integrals of a refractive -index field. In the case of a 
hollow view, as shown in Fig 1, it is feasible to obtain theoretically a 
unique reconstruction of the field outside the opaque object, where the 
information is available.*^ Cha^ and Vest and Prlkryl,^ working on this 
principle, reconstructed incomplete projections by fitting only the available 
data with a truncated series. In practice, however, series expansion methods 
pose a stability problem under severe ill-posed conditions. Past practices 

are mostly to interpolate missing information directly or iteratively with or 
without incorporating a priori information. 

Zien et al.^ assigned fictitious values to the missing zone based on the 

zeroth-moment invariant principle of the line integral transform (Radon 
111 

transform). Lewitt and Bates 11 demonstrated substantial improvements in 
direct reconstruction when missing segments were intemolated with a simole 



curve or a truncated series which satisfied consistency conditions. Ve9t and 
Prikryl^ adopted image-projection-domain revision in their iterative 
convolution reconstruction, which also employed the zeroth -moment invariant 
principle to initally estimate the missing information. This method is in 
essence similar to those of Medoff et al.^ and Choi et al. ^ except that a 
constant field was not assigned to the opaque object region. 

Here we present a different approach for reconstruction of fields with 
opaque objects, which is spirited from the difference impage reconstruction by 
Hefferman and Robb. ^ The method, as described below, is based on iterative 
estimation of the object field and corresponding difference projections only 
in the region where the information is available. 

(1) Make an initial estimation of the object field f Cr , xp > outside the 

SL • 

opaque object. 

(2) Calculate the projection data g (p,0) of the approximate field 

3 

f (r,ij>) only in the region where the measured data g(p,9) is known, 
a 

(3) Find the projections g.(p,9) = g (p,9) - g (p,9) of the difference 

d a 

field f ,(r,tj>) = f(r,i|>) - f where f(r,tj>) is the unknown exact field, 

d a 

(4) Reconstruct the difference field through the inverse Radon 
transformation: f^(r,i|i) = R * {g^(p,0)}. 

(5) Improve the estimate of the object field by f^ + ^(r,tfi) = f*(r,tji) 

3 3 

+ f (r,tji). At this stage, a priori information can be incorporated to further 
d 

refine the new estimate. 

(6) Continue the iteration by returning to step (2) with the new estimate 
unless certain convergence criteria are met. 


The nature of the iterative correction and the use of a priori in- 



formation are similar to those of previous investigations. ^ However, the 
major difference arises in that neither assignment of a constant value in the 
opaque object region nor estimation of missing projections are necessary 
during iteration. The fundamental principle of the algorithm can be perceived 
in a better way by examining the reconstruction error energy. 

In practice, reconstruction from projections does not lead to the 
original field f(x,y) due to imperfect data-aquistion and reconstruction, that 
is, discrete approximate data and finite numerical calculation. Most of these 
operations, however, can be linear or guasi-linear. By assuming a spatially- 
invariant point-spread fuction h(x,y),^ the directly-reconstructed field 
f r (x,y ) is 

f (x,y) ■ h(x,y )** f(x,y) (l) 
r 

where ** denotes two-dimensional convolution. The error engery E is then, by 
applying Rayleigh's theorem,^ 


/*/®|f(x,y) - 


f r (x,y) I 


'dxdy 


= /"/** | [ 1-H(u,v) ]F(u,v) ^du dv 


( 2 ) 


where the upper case denotes the Fourier transform. In the iterative 
algorithm, the difference field is reconstructed. Consequently, the error 



3y comparing Eqs. (2) and (3), we can see that reasonable estimation of the- 
object field can substantially reduce the error energy in the difference field 
method. Fundamentally the proposed algorithm is based on the successive 
estimation of the object field. The application of a priori information in 
step (5) can further reduce the error energy as demonstrated by - Vest and 
Prikryl. 7 

The iterative algorithm has been tested through computer simulation of 
experiments. The initial estimation of the field (step (1)) was made through 
direct reconstruction of the field. For reconstruction, the convolution method 
with the Shepp-Logan filter^ was adopted. The convolution method requires 
full projection data. Linear interpolation was employed for the completion 
of data. Only the boundary constraint, that is, f(x,y) = 0 for |r| > 1, 
was utilized as a priori information. 

Figure 2 in the plot of a four-hump test field given below. 
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This object field is continuous and defined only in the domain of r > 0.4, 
having an opaque object of radius 0.4 at the center. 32 equally-spaced data 
points were generated, excluding the opaque region, for each of 2S projections 

O 

which scanned Che 180 viewing anele# The reconstruction error was defined by 
jff(x,y) - £ (x,y)]/ '-lax { ?(x,v) j J. As shown in Fig. 3, the maximum and 



average errors were reduced from che initial values of 20.9 and 3.6 percent to 
the limiting values of 13.9 and 2.1 percent, respectively, after 12 
iterations. The errors were calculated at 60 x 60 equally-spaced grid points 
in the defined area. The convergence of the errors wa9 relatively slow. 
Figure 4 is the plot of the reconstructed field after 12 iterations. Also 
shown in Fig. 3 are the reconstruction results of the conventional method by 
Vest and Prikryl,^ which iteratively estimates missing projection data through 
projection-object field corrections. The error reduction of the conventional 
method is initially faster than chat of the difference field method; however, 
it diverges. This phenomenon was also observed In the previous 
investigation. ^ 

Reconstruction results from other continuous fields with a simple shape 
also demonstrated a steady but rather slow convergence behavior. The 
convergence is believed to be expedited by employing an appropriate relaxation 
parameter in finding a new approximate object field. The relaxation parameter 
employed in this study is unity. 

In conclusion, this study presented a limited reconstruction example. 
More extensive tests, which cover a broad scope of object fields and 
reconstruction conditions, are due in order to clearly define the behavior of 
the method. As seen in Eq. (3), the convergence of the method depends on the 
initial estimate of th field. One possible way for choosing an appropriate 
estimate is to use the reconstruction from other methods. 
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Figure Captions 


Fig. 1 

Fig. 2 
Fig. 3 
Fig. 4 


Projection of an object field with an opaque object. 
Profile of the test object field. 

Reconstruction errors for two different iteration methods. 
Profile of the reconstructed field after 12 iterations. 












